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Abstract. We investigate the angular and spatial distributions of protons of the 
energy of 0.223 MeV after channeling through an (11, 9) single-wall carbon nanotube 
of the length of 0.2 /xm. The proton incident angle is varied between and 10 mrad, 
being close to the critical angle for channeling. We show that, as the proton incident 
angle increases and approaches the critical angle for channeling, a ring-like structure 
is developed in the angular distribution - donut effect. We demonstrate that it is 
the rainbow effect. When the proton incident angle is between zero and a half of 
the critical angle for channeling, the image force affects considerably the number and 
positions of the maxima of the angular and spatial distributions. However, when the 
proton incident angle is close to the critical angle for channeling, its influence on the 
angular and spatial distributions is reduced strongly. We demonstrate that the increase 
of the proton incident angle can lead to a significant rearrangement of the propagating 
protons within the nanotube. This effect may be used to locate atomic impurities 
in nanotubes as well as for creating nanosized proton beams to be used in materials 
science, biology and medicine. 
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1. Introduction 

While the progress in theoretical modeling and computer simulation of ion channeling 
through carbon nanotubes has reached a mature level, as reviewed in Refs. [T]-[TD], the 
experimental advancement in this area is still in its infancy. Since the issues of ordering, 
straightening and holding nanotubes are probably the most challenging tasks in the 
experimental realization of ion channeling through them, it is not a surprise that the 
best results in performing these tasks are expected when they are grown in a dielectric 
medium. For example, the first experimental data on ion channeling through nanotubes, 
which were reported by Zhu et al. [TT] . were obtained with the He + ions and array of 
well-ordered multi-wall nanotubes grown in a porous anodic aluminum oxide (AI2O3) 
membrane. The authors performed and compared the results of direct measurements 
of the yield of ions transmitted through the bare AI2O3 sample and the AI2O3 sample 
with nanotubes. 

On the other hand, the first experimental results on electron channeling through 
carbon nanotubes were reported by Chai et al. [12]. The authors studied the transport 
of electrons of the energy of 300 keV through the aligned multi-wall nanotubes of the 
lengths of 0.7-3.0 /im embedded in the carbon fiber coatings. The misalignment of 
the nanotubes was up to 1°. Besides, Berdinsky et al. [13J succeeded in growing the 
single-wall carbon nanotubes (SWCNTs) in the ion tracks etched in the S1O2 layers on 
a Si substrate, offering an interesting possibility for the experimental realization of ion 
channeling through nanotubes in a wide range of ion energies. 

Regarding the theoretical modeling and computer simulation of ion channeling 
through carbon nanotubes, we note that the effect of dynamic polarization of the 
nanotube atoms valence electrons by the ion is not usually taken into account pQ- 
[Ti] since its influence at very low and very high energies, of the orders of 1 keV and 
1 GeV, respectively, is negligible. However, it is expected that at medium energies, of 
the order of 1 MeV, this effect contributes significantly to the ion energy loss and gives 
rise to an additional force acting on the ions, called the image force [T5], [16], as it has 
been demonstrated in the computer simulation of the angular distributions of protons 
channeled through the SWCNTs in vacuum [17J. The importance of the image force 
has also been emphasized in the related area of ion transmission through cylindrical 
channels in metals [18] [23] and on ions and molecules moving over supported graphene 
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When the ion channeling dynamics at very low and very high energies is concerned, 
the material surrounding the carbon nanotubes serves predominantly as their passive 
container. However, the ions moving at medium energies induce the strong dynamic 
polarization of both the nanotube atoms valence electrons and the surrounding material, 
which in turn gives rise to a sizeable image force [261 121] • I n these two studies, the image 
force was calculated by a two-dimensional (2D) hydrodynamic model of the nanotube 
atoms valence electrons while the surrounding material was described by a frequency 
dependent dielectric function. On the other hand, the image force has recently been 
shown to influence significantly the rainbow effect in proton channeling through the 
short SWCNTs [2H] and double-wall nanotubes in vacuum [29] as well as through the 
short SWCNTs in the dielectric media pEH Si] - We think that it is important to improve 
our understanding of the role of the image force in the rainbow effect with nanotubes 
because, in analogy with the case of surface ion channeling [32] [31], the measurements of 
this effect can give precise information on both the atomic configuration and interaction 
potentials within nanotubes, which have not yet been explored completely. 

However, in ion channeling experiments, the always present questions are the ones 
of ion beam divergence and misalignment. So, it is important to study the influence of 
the effect of dynamic polarization of carbon nanotubes when the initial ion velocity is 
not parallel to the nanotube axis. Therefore, in this paper, we continue our investigation 
of the image force with the case in which the ion incident angle is not zero. Specifically, 
we analyze the angular and spatial distributions of protons of the velocity of 3 a.u. 
channeled through the straight (11, 9) SWCNTs of the length of 0.2 /im in vacuum. 
The proton incident angle is varied between and 10 mrad, being close to the critical 
angle for channeling. This proton velocity is chosen because the dynamic polarization 
effect is the strongest in the range about it. The consideration is limited to the case 
of a nanotube in vacuum because the presence of a dielectric medium around it would 
introduce only a slight modifying factor in the results of calculation [301 131] . 

It is well known that, for the ion incident angles close to the critical angle for 
channeling, the donut effect develops in the angular distributions of channeled ions. The 
effect was measured with the Si and Ge crystals [35] [37], and explained independently 
afterwards by the theory of crystal rainbows [381 [39] . That theory was formulated 
as a proper theory of ion channeling through thin crystals [40J, and has been applied 
subsequently to ion channeling through short carbon nanotubes [H]-[ll]. It must be 
noted that the donut effect has also been observed in a computer simulation of ion 
propagation through nanotubes [45]. However, the authors did not connect the obtained 
results to the rainbow effect. We explore here the donut effect in the angular and spatial 
distributions of protons channeled through a (11, 9) SWCNT in the presence of the image 
force. 

Regarding the angular and spatial distributions of channeled protons to be 
presented in this study, corresponding to the case in which the proton incident angle is 
not zero, we note that the proton equations of motion in the transverse position plane 



CONTENTS 



4 



that are solved to generate them are 2D. This means that the case we explore is truly 2D, 
unlike the cases treated in our previous studies of the image force in carbon nanotubes, 
which were in fact one- dimensional (ID) [28]~|31j. 

The atomic units will be used throughout the paper unless explicitly stated 
otherwise. 



2. Theory 

We adopt the right Cartesian coordinate system with the z axis coinciding with the 
nanotube axis, the origin in the entrance plane of the nanotube, and the x and y axes 
the vertical and horizontal axes, respectively. The initial proton velocity, v, is taken to 
lie in the yz plane and make angle (p with the z axis, being the proton incident angle. 
The length of the nanotube, L, is assumed to be large enough to allow us to ignore the 
influence of the nanotube edges on the image force, and, at the same time, small enough 
to neglect the energy losses of channeled protons. 

We assume that the interaction between the proton and nanotube atoms can be 
treated classically using the Doyle- Turner expression [16] averaged axially [17] and 
azimuthally [45]. This interaction is repulsive and of the short-range character. Thus, 
the repulsive interaction potential in the proton channeling through the nanotube is of 
the form 



327rdZiZ 2 a 
rcp(r) = 3v/3/ 2 

4 

xJ2 a$I {2ab]r) exp{-6 2 [r 2 + a 2 ]}, (1) 
j=i 

where Z\ — 1 and Z 2 = Q are the atomic numbers of the hydrogen and carbon 
atoms, respectively, a is the nanotube radius, I is the nanotube atoms bond length, 
r = (x 2 + y 2 ) 2 is the distance between the proton and nanotube axis, J is the modified 
Bessel function of the I s * kind and th order, and a, = {0.115,0.188,0.072,0.020} and 
fej = {0.547,0.989, 1.982,5.656} are the fitting parameters (in atomic units) [46J. 

The dynamic polarization of the nanotube by the proton is treated via a 2D 
hydrodynamic model of the nanotube atoms valence electrons, based on a jellium-like 
description of the ion cores making the nanotube wall |15j— 126]. This model includes 
the axial and azimuthal averaging similar to that applied in obtaining the corresponding 
repulsive interaction potential, given by Eq. ([1]). It finally gives the interaction potential 
between the proton and its image, Ui m (r, t), which is stationary in the coordinate system 
moving with the proton and depends on its velocity. This interaction is attractive and 
of the long-range character. The details of derivation of the expression for U- im (r,t) are 
given elsewhere [15] [30]. Consequently, the total interaction potential in the proton 
channeling through the nanotube is 

U(r,t) = U tep (r) + U im (r,t). (2) 
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The proton equations of motion we solve are 
dU(r, t) 

mx ^ = — dx~ ' O 

= — ( 4 ) 

where m is the proton mass. They are subject to the initial conditions for the transverse 
components of the proton velocity that are 

x(t = 0) = 0, (5) 

y(t = 0) = v sin <p « i><£>. (6) 

The longitudinal proton motion is treated as uniform with the initial condition for 
the longitudinal component of the proton velocity that is z(t = 0) = vcostp ~ v. As a 
result, the longitudinal component of the proton position is z(t) = vt. Equations (5) and 
(6) are solved numerically. The angular and spatial distributions of transmitted protons 
are generated using a Monte Carlo computer simulation code. The components of the 
proton impact parameter, xq and yo, are chosen randomly from a uniform distribution 
within the cross-sectional area of the nanotube and its entrance plane. With I = 0.144 
nm |48j, we obtain that a = 0.689 nm. If the proton impact parameter falls inside 
annular interval [a — a sc , a], where a sc = [9n 2 / (128Z2)]3ao is the screening radius and 
the Bohr radius, the proton is treated as if it is backscattered and is disregarded. The 
initial number of protons is about 1 000 000. 

The components of the proton scattering angle, Q x and Q y , are obtained via 
expressions Q x = V x /v and Q y = V y /v, where V x and V y are the final transverse 
components of the proton velocity, which are obtained, together with the final transverse 
components of the proton position, X and Y, as the solutions of Eqs. (jSJ) and (jSJ). The 
proton channeling through the nanotube can be analyzed via the mapping of the impact 
parameter plane, the x y plane, to the scattering angle plane, the Q x Q y plane [40J. 
The corresponding total interaction potential, given by Eq. (j2j), is axially symmetric. 
This means that, if the initial proton velocities were parallel to the nanotube axis, this 
mapping would be ID. However, the initial proton velocities are not parallel to the 
nanotube axis, and the mapping is 2D. Since the proton scattering angle is small, its 
differential transmission cross section is given by 

where Jo is the Jacobian of the mapping, 
_ dQ x dQ y dQ x dQ y 

j0 -&^^-^v (8) 

Thus, equation J® = determines the lines in the impact parameter plane along 
which the proton differential transmission cross section is singular. The images of these 
lines in the scattering angle plane are the rainbow lines in this plane 
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Figure 1. The angular distribution of protons channeled in the (11, 9) SWCNT with 
the inclusion of the image force when the proton incident angle ip = 0. The proton 
velocity is v — 3 a.u., the nanotube radius a — 0.689 nm, and the nanotube length L 
= 0.2 fan. 



We can analyze in a similar way the mapping of the impact parameter plane (the 
xoVo plane), which is the entrance plane of the nanotube and the initial transverse 
position plane, to the exit plane of the nanotube or the final transverse position plane, 
the plane. The Jacobian of this mapping is 

OX dY dXdY 

Jr ~ d^ Wo " Wod^o- (9) 

The rainbow lines in the final transverse position plane are the images of the lines 
in the impact parameter plane determined by equation Jr = 0. 

3. Results and discussion 

Let us now analyze the angular and spatial distributions of protons channeled in the 
(11, 9) SWCNT of the length of 0.2 //m. In all the cases to be studied, the initial 
proton velocity will be v = 3 a.u., corresponding to the initial proton energy of 0.223 
MeV, while the incident proton angle, <p, will be varied between and 10 mrad. The 
maximal proton incident angle will be close to the critical angle for channeling, if> c , being 
about 11 mrad. The analysis will also include the typical proton trajectories through 
the nanotube in the proton phase space. 

In Figs. [TJ|6] we shall display the evolution of the angular distribution of channeled 
protons with the increase of ip. In particular, we shall analyze the development of a 
ring-like structure in the angular distribution under the influence of the image force. 

The scatter plot shown in Fig. [T] represents the angular distribution of channeled 
protons for <p = with the image force included. The corresponding angular distribution 
without the inclusion of the image force contains in its central part only a maximum 
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Figure 2. The angular distribution of protons channeled in the (11, 9) SWCNT with 
the inclusion of the image force when the proton incident angle ip = 6 mrad. The 
proton velocity is v — 3 a.u., the nanotube radius a — 0.689 nm, and the nanotube 
length L = 0.2 /im. 
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Figure 3. The distribution along the <d y axis of protons channeled in the (11, 9) 
SWCNT with the image force taken into account - the solid curve - and without it - 
the dashed curve - when the proton incident angle <p = 6 mrad. The proton velocity 
is v — 3 a.u., the nanotube radius a — 0.689 nm, and the nanotube length L = 0.2 
/im. The former curve corresponds to the angular distribution shown in Fig. 2. 



at the origin. This means that the non-monotonic character of the central part of 
the angular distribution with the image force included is due to the effect of dynamic 
polarization. This was discussed in one of our previous papers [28]. We presented in it 
the distributions of channeled protons along the Q y axis (in the scattering angle plane) 
with and without the image force included, which were in fact (for if = 0) the radial 
yields of channeled protons. The conclusion of the discussion was that the maxima of 
the radial yield, appearing when the image force was included, were due to the rainbow 
effect. 
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Figure 4. The angular distribution of protons channeled in the (11, 9) SWCNT with 
the inclusion of the image force when the proton incident angle tp = 10 mrad. The 
proton velocity is v — 3 a.u., the nanotube radius a — 0.689 nm, and the nanotube 
length L = 0.2 fim. 
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Figure 5. 2The distribution along the Oj, axis of protons channeled in the (11, 9) 
SWCNT with the image force taken into account - the solid curve - and without it - 
the dashed curve - when the proton incident angle tp — 10 mrad. The proton velocity 
is v = 3 a.u., the nanotube radius a = 0.689 nm, and the nanotube length L = 0.2 
/im. The former curve corresponds to the angular distribution shown in Fig. 4. 

Fig. |2] gives the angular distribution of channeled protons for ip = 6 mrad with the 
image force included. One can notice easily about a half of a ring-like structure, with an 
exceptionally high yield of channeled protons. This is the precursor of the effect known 
as the donut effect, which is connected to the misalignment of the proton beam and 
nanotube axis. In addition, the angular distribution contains several intricately shaped 
regions with lower yields of channeled protons. We show in Fig. |3] the corresponding 
distribution of channeled protons along the Q y axis with and without the image force 
included. The sharp maximum of this distribution, appearing at -6.0 mrad, is due to 
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Figure 6. The rainbow lines in the scattering angle plane for the protons channeled 
in the (11, 9) SWCNT with the inclusion of the image force when the proton incident 
angle tp = 10 mrad. The proton velocity is v — 3 a.u., the nanotube radius a — 0.689 
nm, and the nanotube length L = 0.2 fim. 
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Figure 7. The spatial distribution of protons channeled in the (11, 9) SWCNT with 
the image force included when the proton incident angle tp = 0. The proton velocity is 
v — 3 a.u., the nanotube radius a — 0.689 nm, and the nanotube length L = 0.2 fjxa. 

the donut effect. It is evident that the image force makes this maximum weaker. On 
the other hand, the origin of the broad maximum of the distribution, located at -5.1 
mrad, is solely the image force. Thus, we can conclude that for the median values of 
tp, between and about ip c /2, the image force still plays a significant role in generating 
the angular distribution. 

We show in Fig. HJthe angular distribution of channeled protons for tp = 10 mrad 
with the image force included. One can see clearly the whole ring-like structure, with 
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Figure 8. The spatial distribution of protons channeled in the (11, 9) SWCNT with 
the inclusion of the image force when the proton incident angle tp — 10 mrad. The 
proton velocity is v — 3 a.u., the nanotube radius a — 0.689 nm, and the nanotube 
length L = 0.2 [im. 
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Figure 9. The distribution along the Y axis of protons channeled in the (11, 9) 
SWCNT with the image force taken into account - the solid curve - and without it - 
the dashed curve - when the proton incident angle (p = 10 mrad. The proton velocity 
is v = 3 a.u., the nanotube radius a = 0.689 nm, and the nanotube length L = 0.2 
(im. The former curve corresponds to the spatial distribution given in Fig. 8. 

an exceptionally high yield of channeled protons. This is the fully developed donut 
effect. As it has been already said, the corresponding value of (p is close to the value 
of ip c . In addition, as in Fig. [21 the angular distribution contains several intricately 
shaped regions with lower but distinctly graded yields of channeled protons, with the 
very clear boundaries between them. Fig. [5] gives the corresponding distribution of 
channeled protons along the Q y axis with and without the image force included. It is 
evident that, when (p is close to ip c , the role of the effect of dynamic polarization in 
generating the angular distribution is almost negligible. Fig. [6] shows the corresponding 
rainbow lines in the scattering angle plane with the dynamic polarization effect taken 
into account. These lines clearly demonstrate that the non- uniformity of the angular 
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Figure 10. The rainbow lines in the final transverse position plane for the protons 
channeled in the (11, 9) SWCNT with the image force included when the proton 
incident angle ip = 10 mrad. The proton velocity is v = 3 a.u., the nanotube radius a 
= 0.689 nm, and the nanotube length L = 0.2 /um. 
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Figure 11. The dependence of the y component of the proton scattering angle on the 
z component of its position in the channeling through the (11, 9) SWCNT with the 
inclusion of the image force when the proton incident angle tp = for the x component 
of the proton impact parameter x = and the y component of its impact parameter 
ya = ±2, ±6 and ±10 a.u. The proton velocity is v = 3 a.u., the nanotube radius a = 
0.689 nm, and the nanotube length L = 0.2 /im. 

distribution, including the donut effect, is due to the rainbow effect. 

In Figs. ITlfTUI we shall display the evolution of the spatial distribution of channeled 
protons with the increase of tp, which is going on in parallel with the evolution of the 
angular distribution displayed in Figs. dH 

The scatter plot given in Fig. [7] represents the spatial distribution of channeled 
protons for ip = with the image force included. This spatial distribution and the 
corresponding spatial distribution without the image force included were analyzed in 
one of our previous papers [51] . We presented in it the distributions of channeled 
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Figure 12. The dependence of the y component of the proton position on the z 
component of its position in the channeling through the (11, 9) SWCNT with the 
inclusion of the image force when the proton incident angle ip = for the x component 
of the proton impact parameter x n = and the y component of its impact parameter 
yo = ±2, ±6 and ±10 a.u. The proton velocity is v — 3 a.u., the nanotube radius = 
0.689 nm, and the nanotube length L = 0.2 /im. 
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Figure 13. The dependence of the y component of the proton scattering angle on the 
z component of its position in the channeling through the (11, 9) SWCNT with the 
inclusion of the image force when the proton incident angle ip — 10 mrad for the x 
component of the proton impact parameter xq = and the y component of its impact 
parameter y = ±2, ±6 and ±10 a.u. The proton velocity is v — 3 a.u., the nanotube 
radius a = 0.689 nm, and the nanotube length L = 0.2 /jm. 



protons along the y axis with and without the image force included, which were in fact 
(for ip = 0) the radial yields of channeled protons. It was demonstrated that the maxima 
of the radial yields, present in both spatial distributions, were the rainbow maxima. We 
also concluded that the dynamic polarization effect caused the shifts of the maxima of 
the spatial distribution generated with the effect not taken into account as well as the 
appearance of the additional maxima. 

The spatial distribution of channeled protons for ip — 10 mrad with the effect of 
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Figure 14. The dependence of the y component of the proton position on the z 
component of its position in the channeling through the (11, 9) SWCNT with the 
image force taken into account when the proton incident angle ip = 10 mrad for the x 
component of the proton impact parameter xq — and the y component of its impact 
parameter yo = ±2, ±6 and ±10 a.u. The proton velocity is v — 3 a.u., the nanotube 
radius a = 0.689 nm, and the nanotube length L = 0.2 fim. 

dynamic polarization taken into account is presented in Fig. [HJ When this spatial 
distribution is compared to the spatial distribution for (p — 0, it is evident that the 
maximal change of ip induces a significant rearrangement of the protons in the final 
transverse position plane. Almost all the protons are displaced to the left half of the 
nanotube. However, as in the cases of angular distributions for (p — 6 and 10 mrad, 
the spatial distribution also contains several intricately shaped regions with lower but 
distinctly graded yields of channeled protons, with the very clear boundaries between 
them. Fig. |9] gives the corresponding distribution of channeled protons along the Y 
axis, in the final transverse position plane, with and without the image force included. 
For this value of ip, the strongest maximum of the spatial distribution lies at -3.6 a.u. 
instead at the origin for ip = 0. One can also conclude that, when ip is close to ip c , the 
role of the effect of dynamic polarization in generating the spatial distribution is small 
but noticeable. The effect makes the strongest maximum of the spatial distribution 
weaker and induces a rightward shift of the second strongest maximum. Fig. [10] shows 
the corresponding rainbow lines in the final transverse position plane with the image 
force taken into account. As in the case of angular distribution for this value of (p, these 
lines clearly demonstrate that the non-uniformity of the spatial distribution is to be 
attributed to the rainbow effect. 

In Figs. [Tlin^l we shall display the typical proton trajectories through the nanotube 
in the proton phase space, complementing the result displayed in Figs. ITHTOl 

We show in Fig. [11] the y component of the proton scattering angle (0^) as a 
function of the z component of its position within the nanotube with the effect of 
dynamic polarization included when p = for xq = and y = ±2, ±6 and ±10 a.u. 
Looking at the angular distribution shown in Fig. [T] we see that the channeled protons 
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with the impact parameters close to the nanotube axis, i.e., for y = ±2 a.u., and to the 
nanotube wall, i.e., for y = ±10 a.u., contribute to the part of the angular distribution 
close to the origin. The channeled protons with the impact parameters comparable to 
a/2, i.e., for y = ±6 a.u., give rise to the rainbow maxima lying at about 2 mrad. 

Fig. [12] gives the dependence of the y component of the proton position on the z 
component of its position within the nanotube with the image force included when (p = 
for the same values of the components of the proton impact parameter as in Fig. [TTJ 
One can see that the channeled protons with yo = ±2 a.u. give rise to the part of the 
spatial distribution close to the origin. The channeled protons with y = ±6 and ±10 
a.u. contribute to the peripheral part of the spatial distribution. 

We give in Fig. [13]the y component of the proton scattering angle (0^) as a function 
of the z component of its position with the image force taken into account when <p = 
10 mrad for the same values of the components of the proton impact parameter as in 
Fig. [HJ It is easy to conclude that the channeled protons with y Q = 2, 6 and 10 a.u. 
contribute to the right part of the donut. The channeled protons with yo = -2, -6 and 
-10 a.u. give rise to the most intense part of the donut, being its farthest left part. 

Fig. E] shows the dependence of the y component of the proton position on the z 
component of its position with the image force included when ip = 10 mrad for the same 
values of the components of the proton impact parameter as in Fig. [TTJ It is evident 
that all the propagating protons in question end up in the left half of the nanotube, 
after being reflected from the right part of the nanotube wall. 

An additional result of our computer simulations is related to the influence of the 
image force on tp c . We followed the change of the total yield of channeled protons with 
the increase of <p>, and found that with the image force taken into account ip c = 10.6 
mrad. When the image force is not taken into account ip c = 11.9 mrad. This increase 
of "0c is attributed to the increase of the total interaction potential in question, given 
by Eq. (j2J), when its attractive component, originating in the interaction of the proton 
and its image, is not taken into account. This conclusion is justified via the relation 
between %p c and the total interaction potential at the distance from the nanotube wall 
equal to the screening radius, Use, 



where E is the initial proton energy |49j . 
4. Concluding remarks 

We have presented the first theoretical investigation of the angular and spatial 
distributions of ions channeled through a nanotube for different proton incidence angles 
with the effect of dynamic polarization of the nanotube included. The ions are protons 
of the velocity of v = 3 a.u. and the nanotube is an (11, 9) SWCNT of the length of L 
= 0.2 /iin. The proton incident angle, ip, is varied between and 10 mrad, being close 




(10) 
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to the critical angle for channeling, ip c . We have noticed a slight increase of ip c when 
the image force is not taken into account. 

We have observed a ring-like structure developing in the angular distribution of 
channeled protons with tp increasing and approaching ip c . The effect has been recognized 
as the donut effect, being in fact the rainbow effect. If <p is between and about ip c /2, the 
image force plays a significant role in generating the angular and spatial distributions, 
including the rainbow maxima. However, if ip is close to ip c , the contribution of the image 
force to the angular and spatial distributions, including the donut effect, is minor. 

The analysis of the generated spatial distributions of channeled protons has shown 
that the increase of tp can give rise to a significant rearrangement of the propagating 
protons within the nanotube. For example, for tp — 10 mrad, the proton beam is 
displaced from the nanotube axis toward the nanotube wall leaving the region around the 
axis practically empty. It is clear that such a rearrangement of the propagating protons 
may be used to locate various atomic impurities in the nanotube, using the secondary 
processes like backward Coulomb scattering and nuclear reactions. In addition, the 
presence of the rainbow maxima in the spatial distributions can be used to determine 
the positions of the impurities very precisely. One can also think about directing such a 
nonosized proton beam to a material to be modified with it, or to a biological or medical 
sample. 
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